home *** CD-ROM | disk | FTP | other *** search
/ The Atari Compendium / The Atari Compendium (Toad Computers) (1994).iso / files / prgtools / gnustuff / tos / updates / update16.zoo / pml / diffs
Encoding:
Text File  |  1992-03-06  |  21.0 KB  |  992 lines
  1. *** 1.12    1992/02/04 16:42:12
  2. --- Changelog    1992/03/06 19:37:53
  3. ***************
  4. *** 252,254 ****
  5. --- 252,270 ----
  6.   the correct parameters.
  7.   
  8.   ---------------------------- Patchlevel 15 ------------------------------
  9. + 
  10. + Makefile.tos:: Michal
  11. +     use macros to make it easier to to adjust.
  12. + 
  13. + atan2.c: mjr++
  14. + ccos.c:
  15. + cdiv.c:
  16. + cmult.c:
  17. + csin.c:
  18. +     new assembler routines for both the sfp004 and the 881/882
  19. + 
  20. + cadd.c: mjr++
  21. + csubt.c:
  22. +     fixed nasty arg passing bug in the sfp code.
  23. + 
  24. + ---------------------------- Patchlevel 16 ------------------------------
  25. *** 1.7    1992/02/03 20:22:27
  26. --- Makefile.tos    1992/03/06 19:37:54
  27. ***************
  28. *** 3,8 ****
  29. --- 3,12 ----
  30.   CROSSBIN = $(CROSSDIR)/bin
  31.   CROSSINC = $(CROSSDIR)/include
  32.   XFLAGS = -Wstrict-prototypes
  33. + #this has to be gnu make
  34. + MAKE=make
  35. + CC=cgcc
  36. + AR=car
  37.   
  38.   ALL = pml16.olb pml.olb
  39.   
  40. ***************
  41. *** 9,28 ****
  42.   all : $(ALL)
  43.   
  44.   pml16.olb :
  45. !     make -f makefile.pml clean
  46. !     make -f makefile.pml CC="cgcc" AR="car" XFLAGS="$(XFLAGS)" s=true
  47.   
  48.   pml.olb :
  49. !     make -f makefile.pml clean
  50. !     make -f makefile.pml CC="cgcc" AR="car" XFLAGS="$(XFLAGS)"
  51.   
  52.   pml02016.olb :
  53. !     make -f makefile.pml clean
  54. !     make -f makefile.pml CC="cgcc -m68020" AR="car" XFLAGS="$(XFLAGS)" m=68881 s=true
  55.   
  56.   pml020.olb :
  57. !     make -f makefile.pml clean
  58. !     make -f makefile.pml CC="cgcc -m68020" AR="car" XFLAGS="$(XFLAGS)" m=68881
  59.   
  60.   install : $(ALL)
  61.       cp pml16.olb pml.olb $(CROSSLIB)
  62. --- 13,34 ----
  63.   all : $(ALL)
  64.   
  65.   pml16.olb :
  66. !     $(MAKE) -f makefile.pml clean
  67. !     $(MAKE) -f makefile.pml CC="$(CC)" AR="$(AR)" XFLAGS="$(XFLAGS)" s=true
  68.   
  69.   pml.olb :
  70. !     $(MAKE) -f makefile.pml clean
  71. !     $(MAKE) -f makefile.pml CC="$(CC)" AR="$(AR)" XFLAGS="$(XFLAGS)"
  72.   
  73.   pml02016.olb :
  74. !     $(MAKE) -f makefile.pml clean
  75. !     $(MAKE) -f makefile.pml CC="$(CC) -m68020" AR="$(AR)" \
  76. !  XFLAGS="$(XFLAGS)" m=68881 s=true
  77.   
  78.   pml020.olb :
  79. !     $(MAKE) -f makefile.pml clean
  80. !     $(MAKE) -f makefile.pml CC="$(CC) -m68020" AR="$(AR)" \
  81. !  XFLAGS="$(XFLAGS)" m=68881
  82.   
  83.   install : $(ALL)
  84.       cp pml16.olb pml.olb $(CROSSLIB)
  85. ***************
  86. *** 33,41 ****
  87.       cp math.h $(CROSSINC)
  88.   
  89.   clean:
  90. !     make -f makefile.pml clean
  91. !     make -f makefile.pml clean
  92.   
  93.   realclean:
  94. !     make -f makefile.pml realclean
  95. !     make -f makefile.pml realclean
  96. --- 39,47 ----
  97.       cp math.h $(CROSSINC)
  98.   
  99.   clean:
  100. !     $(MAKE) -f makefile.pml clean
  101. !     $(MAKE) -f makefile.pml clean
  102.   
  103.   realclean:
  104. !     $(MAKE) -f makefile.pml realclean
  105. !     $(MAKE) -f makefile.pml realclean
  106. *** 1.12    1992/02/04 16:42:12
  107. --- PatchLevel.h    1992/03/06 19:37:54
  108. ***************
  109. *** 1,5 ****
  110.   /*
  111. !  *    PatchLevel: 15
  112.    *
  113.    *    the Patch Level above is to identify the version
  114.    *    of the all the files in this directory. given the above
  115. --- 1,5 ----
  116.   /*
  117. !  *    PatchLevel: 16
  118.    *
  119.    *    the Patch Level above is to identify the version
  120.    *    of the all the files in this directory. given the above
  121. *** 1.7    1992/02/03 20:19:23
  122. --- atan2.c    1992/03/06 19:37:56
  123. ***************
  124. *** 87,92 ****
  125. --- 87,94 ----
  126.    *
  127.    */
  128.   
  129. + #if !defined (__M68881__) && !defined (sfp004)
  130. + 
  131.   #include <stdio.h>
  132.   #include <math.h>
  133.   #include "pml.h"
  134. ***************
  135. *** 124,126 ****
  136. --- 126,280 ----
  137.       return (result);
  138.   }
  139.   
  140. + #endif !defined (__M68881__) && !defined (sfp004)
  141. + #ifdef    __M68881__
  142. + __asm("
  143. + .text
  144. + .even
  145. + _funcname:
  146. +     .ascii    \"atan2\\0\"
  147. +     .even
  148. + 
  149. + .globl    _atan2
  150. + _atan2:
  151. + | denormalized numbers are treated as 0
  152. +     tstl    sp@(12)
  153. +     beq    5f        | x == 0!
  154. +     blt    1f        | x < 0!
  155. +                 | x > 0: return atan(y/x)
  156. + 
  157. +     fmoved    sp@(4)@,fp0    | get y
  158. +     fdivd    sp@(12),fp0    | y/x    
  159. +     fatanx    fp0,fp0        | atan(y/x)
  160. +     bra 3f            | return
  161. + 1:                | x < 0
  162. + 
  163. +     fmovecr    #0,fp1        | get pi
  164. +     fmoved    sp@(4)@,fp0    | get y
  165. +     fdivd    sp@(12),fp0    | y/x
  166. +     fatanx    fp0,fp0        | atan(y/x)
  167. +     btst    #31,sp@(4)    | sign(y)
  168. +     beq    2f        | positive!
  169. + 
  170. +     fnegx    fp1,fp1        | transfer sign
  171. + 2:    faddx    fp1,fp0        | sign(y)*pi + atan(y/x)
  172. + |    bra 3f            | return
  173. + 3:
  174. +     fmoved    fp0,sp@-    | return result
  175. +     moveml    sp@+,d0/d1
  176. + 4:    
  177. +     rts            | sigh.
  178. + 5:                | x == 0
  179. +     movel    #1073291771,d0    | pi/2
  180. +     movel    #1413754136,d1    |
  181. + 
  182. +     tstl    sp@(4)        | 
  183. +     beq    6f        | NaN
  184. +     bge    4b        | exit
  185. +     bset    #31,d0        | x < 0 : return -pi/2
  186. +     bra    4b
  187. + 6:    movel    #-1,d0        | NaN
  188. +     movel    #-1,d1        |
  189. +     bra    4b
  190. + ");    /* end asm    */
  191. + #endif    __M68881__
  192. + 
  193. + #ifdef    sfp004
  194. + __asm("
  195. + 
  196. + comm =     -6
  197. + resp =    -16
  198. + zahl =      0
  199. + 
  200. + .even
  201. + .text
  202. + .even
  203. + _funcname:
  204. +     .ascii    \"atan2\\0\"
  205. +     .even
  206. + .text
  207. + .even
  208. + .globl    _atan2
  209. + _atan2:
  210. + | denormalized numbers are treated as 0
  211. +     lea    0xfffa50,a0
  212. +     moveml    a7@(12),d0-d1    |  x
  213. +     tstl    d0
  214. +     beq    5f        | x == 0!
  215. +     blt    1f        | x < 0!
  216. +                 | x > 0: return atan(y/x)
  217. + 
  218. + |    fmoved    sp@(4),fp0    | get y
  219. +     movew    #0x5400,a0@(comm)
  220. +     .long    0x0c688900, 0xfff067f8
  221. +     movel    sp@(4),a0@
  222. +     movel    sp@(8),a0@
  223. + 
  224. + |    fdivd    sp@(12),fp0    | y/x
  225. +     movew    #0x5420,a0@(comm)
  226. +     .long    0x0c688900, 0xfff067f8
  227. +     movel    d0,a0@
  228. +     movel    d1,a0@
  229. + 
  230. + |    fatanx    fp0,fp0        | atan(y/x)
  231. +     movew    #0x000a,a0@(comm)
  232. +     .word    0x4a68,0xfff0,0x6bfa
  233. + 
  234. +     bra 3f            | return
  235. + 1:                | x < 0
  236. + 
  237. + |    fmovecr    #0,fp1        | get pi
  238. +     movew    #0x5c80,a0@(comm)
  239. +     .long    0x0c688900, 0xfff067f8
  240. + 
  241. + |    fmoved    sp@(4),fp0    | get y
  242. +     movew    #0x5400,a0@(comm)
  243. +     .long    0x0c688900, 0xfff067f8
  244. +     movel    sp@(4),a0@
  245. +     movel    sp@(8),a0@
  246. + 
  247. + |    fdivd    sp@(12),fp0    | y/x
  248. +     movew    #0x5420,a0@(comm)
  249. +     .long    0x0c688900, 0xfff067f8
  250. +     movel    d0,a0@
  251. +     movel    d1,a0@
  252. + 
  253. + |    fatanx    fp0,fp0        | atan(y/x)
  254. +     movew    #0x000a,a0@(comm)
  255. +     .word    0x4a68,0xfff0,0x6bfa
  256. + 
  257. +     btst    #31,sp@(4)    | sign(y)
  258. +     beq    2f        | positive!
  259. + 
  260. + |    fnegx    fp1,fp1        | transfer sign
  261. +     movew    #0x049a,a0@(comm)
  262. +     .word    0x4a68,0xfff0,0x6bfa
  263. + 
  264. + 2:|    faddx    fp1,fp0        | sign(y)*pi + atan(y/x)
  265. +     movew    #0x0422,a0@(comm)
  266. +     .word    0x4a68,0xfff0,0x6bfa
  267. + 
  268. + |    bra 3f            | return
  269. + 3:
  270. + |    fmoved    fp0,d0-d1    | return result
  271. +     movew    #0x7400,a0@(comm)
  272. +     .long    0x0c688900, 0xfff067f8
  273. +     movel    a0@,d0
  274. +     movel    a0@,d1
  275. + 
  276. + 4:    
  277. +     rts            | sigh.
  278. + 5:                | x == 0
  279. +     movel    #1073291771,d0    | pi/2
  280. +     movel    #1413754136,d1    |
  281. + 
  282. +     tstl    sp@(4)        | 
  283. +     beq    6f        | NaN
  284. +     bge    4b        | exit
  285. +     bset    #31,d0        | x < 0 : return -pi/2
  286. +     bra    4b
  287. + 6:    movel    #-1,d0        | NaN
  288. +     movel    #-1,d1        |
  289. +     bra    4b
  290. + ");    /* end asm    */
  291. + #endif    sfp004
  292. *** 1.3    1992/02/03 20:19:23
  293. --- cadd.c    1992/03/06 19:37:57
  294. ***************
  295. *** 120,127 ****
  296.   
  297.       movew    #0x5422,a0@(comm)    | fp0 += z2.real
  298.       .long    0x0c688900, 0xfff067f8
  299. !     movel    a7@(4),a0@        | load arg_hi
  300. !     movel    a7@(8),a0@        | load arg_low
  301.   
  302.       movew    #0x5480,a0@(comm)    | z1.imag -> fp1
  303.       .long    0x0c688900, 0xfff067f8
  304. --- 120,127 ----
  305.   
  306.       movew    #0x5422,a0@(comm)    | fp0 += z2.real
  307.       .long    0x0c688900, 0xfff067f8
  308. !     movel    a7@(20),a0@        | load arg_hi
  309. !     movel    a7@(24),a0@        | load arg_low
  310.   
  311.       movew    #0x5480,a0@(comm)    | z1.imag -> fp1
  312.       .long    0x0c688900, 0xfff067f8
  313. ***************
  314. *** 129,138 ****
  315.       movel    a7@(16),a0@        | load arg_low
  316.   
  317.       movew    #0x54a2,a0@(comm)    | fp1 += z2.imag
  318. !     movel    a1,d0             pointer to result
  319.       .long    0x0c688900, 0xfff067f8
  320. !     movel    a7@(4),a0@        | load arg_hi
  321. !     movel    a7@(8),a0@        | load arg_low
  322.   
  323.   |    movew    #%0111 0101 0000 0000,a0@(comm)    | fetch fp0
  324.       movew    #0x7400,a0@(comm)        | 
  325. --- 129,138 ----
  326.       movel    a7@(16),a0@        | load arg_low
  327.   
  328.       movew    #0x54a2,a0@(comm)    | fp1 += z2.imag
  329. !     movel    a1,d0            | pointer to result
  330.       .long    0x0c688900, 0xfff067f8
  331. !     movel    a7@(28),a0@        | load arg_hi
  332. !     movel    a7@(32),a0@        | load arg_low
  333.   
  334.   |    movew    #%0111 0101 0000 0000,a0@(comm)    | fetch fp0
  335.       movew    #0x7400,a0@(comm)        | 
  336. *** 1.4    1992/02/03 20:19:23
  337. --- ccos.c    1992/03/06 19:37:58
  338. ***************
  339. *** 61,69 ****
  340.    *
  341.    */
  342.   
  343. ! #if defined (__M68881__) && !defined (_M68881)
  344. ! /*# define _M68881*/
  345. ! #endif
  346.   
  347.   #include <stdio.h>
  348.   #include <math.h>
  349. --- 61,67 ----
  350.    *
  351.    */
  352.   
  353. ! #if !defined (__M68881__) && !defined (sfp004)
  354.   
  355.   #include <stdio.h>
  356.   #include <math.h>
  357. ***************
  358. *** 79,81 ****
  359. --- 77,195 ----
  360.       result.imag = -sin(z.real) * sinh(z.imag);
  361.       return (result);
  362.   }
  363. + #endif !defined (__M68881__) && !defined (sfp004)
  364. + 
  365. + #ifdef    __M68881__
  366. + __asm("
  367. + .text
  368. + _funcname:
  369. +     .ascii    \"ccos\\0\"
  370. +     .even
  371. +     .globl _ccos
  372. + _ccos:
  373. +     movel    a1,d0        | save a1 as return value
  374. +     fmoved    sp@(4),fp0    | z.real
  375. +     fsinx    fp0,fp1        | cos(z.real)
  376. +     fcosx    fp0,fp0        | sin(z.real)
  377. +     fnegx    fp1,fp1        |-sin(z.real)
  378. + 
  379. +     fmoved    sp@(12),fp3    | z.imag
  380. +     fcoshx    fp3,fp2        | cosh(z.imag)
  381. +     fsinhx    fp3,fp3        | sinh(z.imag)
  382. +     
  383. +     fmulx    fp2,fp0        | result.real
  384. +     fmoved    fp0,a1@        |
  385. +     fmulx    fp3,fp1        | result.imag
  386. +     fmoved    fp1,a1@(8)    |
  387. + ");    /* end asm    */
  388. + #endif    __M68881__
  389. + 
  390. + #ifdef    sfp004
  391. + __asm("
  392. + | double precision floating point stuff for Atari-gcc using the SFP004
  393. + | developed with gas
  394. + |
  395. + | double precision complex sin
  396. + |
  397. + | M. Ritzert (mjr at dmzrzu71)
  398. + |
  399. + | 12.10.1990
  400. + |
  401. + | addresses of the 68881 data port. This choice is fastest when much data is
  402. + | transferred between the two processors.
  403. + 
  404. + comm =     -6
  405. + resp =    -16
  406. + zahl =      0
  407. + 
  408. + | waiting loop ...
  409. + |
  410. + | wait:
  411. + | ww:    cmpiw    #0x8900,a1@(resp)
  412. + |     beq    ww
  413. + | is coded directly by
  414. + |    .long    0x0c688900, 0xfff067f8
  415. + | and
  416. + | www:    tst.b    a1@(resp)
  417. + |    bmi.b    www
  418. + | is coded by
  419. + |    .word    0x4a68,0xfff0,0x6bfa        | test
  420. + 
  421. +     .text; .even
  422. +     .globl _ccos
  423. + _ccos:
  424. +     movel    a1,d0                | save a1 as return value
  425. +     lea    0xfffa50,a0            | fpu address
  426. + 
  427. +     movew    #0x54b0,a0@(comm)        | sincos: sin -> fp1
  428. +     .long    0x0c688900, 0xfff067f8        |      cos -> fp0
  429. +     movel    sp@(4), a0@            | load z.real
  430. +     movel    sp@(8), a0@            | load z.real
  431. + 
  432. + |    fmoved    sp@(12),fp3            | z.imag to fp3
  433. +     movew    #0x5580,a0@(comm)
  434. +     .long    0x0c688900, 0xfff067f8
  435. +     movel    sp@(12),a0@
  436. +     movel    sp@(16),a0@
  437. + 
  438. + |    fcoshx    fp3,fp2                | cosh(z.imag)
  439. +     movew    #0x0d19,a0@(comm)
  440. +     .word    0x4a68,0xfff0,0x6bfa        | test
  441. + |    fsinhx    fp3,fp3                | sinh(z.imag)
  442. +     movew    #0x0d82,a0@(comm)
  443. +     .word    0x4a68,0xfff0,0x6bfa        | test
  444. + |    fmulx    fp2,fp0                | result.real
  445. +     movew    #0x0823,a0@(comm)
  446. +     .word    0x4a68,0xfff0,0x6bfa        | test
  447. + |    fmoved    fp0,a1@                |
  448. +     movew    #0x7400,a0@(comm)        | 
  449. +     .long    0x0c688900, 0xfff067f8
  450. +     movel    a0@,a1@
  451. +     movel    a0@,a1@(4)
  452. + |    fmulx    fp3,fp1                | result.imag
  453. +     movew    #0x0ca3,a0@(comm)
  454. +     .word    0x4a68,0xfff0,0x6bfa        | test
  455. + |    fnegx    fp1,fp1
  456. +     movew    #0x049a,a0@(comm)
  457. +     .word    0x4a68,0xfff0,0x6bfa        | test
  458. + |    fmoves    fp1,d1                |
  459. +     movew    #0x7480,a0@(comm)        |
  460. +     .long    0x0c688900, 0xfff067f8
  461. +     movel    a0@,a1@(8)
  462. +     movel    a0@,a1@(12)
  463. + ");    /* end asm    */
  464. + #endif    sfp004
  465. + 
  466. + 
  467. + #if defined (__M68881__) || defined (sfp004)
  468. + # ifdef ERROR_CHECK    /* no error checking for now    */
  469. + __asm("
  470. +     pea    _funcname
  471. +     jmp    c_err_check
  472. + ");    /* end asm    */
  473. + # else  ERROR_CHECK
  474. + 
  475. + __asm("rts");
  476. + 
  477. + # endif ERROR_CHECK
  478. + #endif defined (__M68881__) || defined (sfp004)
  479. *** 1.3    1992/02/03 20:19:23
  480. --- cdiv.c    1992/03/06 19:37:59
  481. ***************
  482. *** 68,76 ****
  483.    *
  484.    */
  485.   
  486. ! #if defined (__M68881__) && !defined (_M68881)
  487. ! /*# define _M68881*/
  488. ! #endif
  489.   
  490.   #include <stdio.h>
  491.   #include <math.h>
  492. --- 68,74 ----
  493.    *
  494.    */
  495.   
  496. ! #if !defined (__M68881__) && !defined (sfp004)
  497.   
  498.   #include <stdio.h>
  499.   #include <math.h>
  500. ***************
  501. *** 114,116 ****
  502. --- 112,268 ----
  503.       }
  504.       return (result);
  505.   }
  506. + #endif !defined (__M68881__) && !defined (sfp004)
  507. + #ifdef    __M68881__
  508. + __asm("
  509. + .text
  510. + .even
  511. + _funcname:
  512. +     .ascii    \"cdiv\\0\"
  513. +     .even
  514. + 
  515. + .globl    _cdiv
  516. + _cdiv:
  517. +     fmoved    sp@(4),fp0
  518. +     fmoved    sp@(12),fp1
  519. +     fmoved    sp@(20),fp2
  520. +     fmoved    sp@(28),fp3
  521. +     fmovex    fp0,fp4
  522. +     movel    a1,d0        | pointer to result
  523. + 
  524. +     fmovex    fp2,fp5
  525. +     fmulx    fp2,fp5
  526. +     fmovex    fp3,fp6
  527. +     fmulx    fp3,fp6
  528. +     faddx    fp6,fp5
  529. + 
  530. +     fmulx    fp2,fp4
  531. +     fmulx    fp3,fp0
  532. +     fmulx    fp1,fp2    
  533. +     fmulx    fp1,fp3
  534. +     faddx    fp3,fp4
  535. +     fdivx    fp5,fp4
  536. +     fsubx    fp0,fp2
  537. +     fdivx    fp5,fp2
  538. + 
  539. +     fmoved    fp4,a1@
  540. +     fmoved    fp2,a1@(8)
  541. + ");    /* end asm    */
  542. + #endif    __M68881__
  543. + 
  544. + #ifdef    sfp004
  545. + __asm("
  546. + 
  547. + comm =     -6
  548. + resp =    -16
  549. + zahl =      0
  550. + 
  551. + .even
  552. + .text
  553. + .even
  554. + _funcname:
  555. +     .ascii    \"cdiv\\0\"
  556. +     .even
  557. + .text
  558. + .even
  559. + .globl    _cdiv
  560. + _cdiv:
  561. + 
  562. +     lea    0xfffa50,a0
  563. + 
  564. +     movew    #0x5400,a0@(comm)    | z1.real -> fp0
  565. +     movel    a1,d0            | pointer to result
  566. +     .long    0x0c688900, 0xfff067f8
  567. +     movel    a7@(4),a0@        | load arg_hi
  568. +     movel    a7@(8),a0@        | load arg_low
  569. + 
  570. +     movew    #0x5480,a0@(comm)    | z1.imag -> fp1
  571. +     .long    0x0c688900, 0xfff067f8
  572. +     movel    a7@(12),a0@        | load arg_hi
  573. +     movel    a7@(16),a0@        | load arg_low
  574. + 
  575. +     movew    #0x5500,a0@(comm)    | z2.real -> fp2
  576. +     .long    0x0c688900, 0xfff067f8
  577. +     movel    a7@(20),a0@        | load arg_hi
  578. +     movel    a7@(24),a0@        | load arg_low
  579. + 
  580. +     movew    #0x5580,a0@(comm)    | z2.imag -> fp3
  581. +     .long    0x0c688900, 0xfff067f8
  582. +     movel    a7@(28),a0@        | load arg_hi
  583. +     movel    a7@(32),a0@        | load arg_low
  584. + 
  585. +     movew    #0x0200,a0@(comm)    | copy fp0 to fp4
  586. +     .word    0x4a68,0xfff0,0x6bfa    | test
  587. + 
  588. + 
  589. + |    fmovex    fp2,fp5
  590. +     movew    #0x0a80,a0@(comm)
  591. +     .word    0x4a68,0xfff0,0x6bfa    | test
  592. + |    fmulx    fp2,fp5
  593. +     movew    #0x0aa3,a0@(comm)
  594. +     .word    0x4a68,0xfff0,0x6bfa    | test
  595. + |    fmovex    fp3,fp6
  596. +     movew    #0x0f00,a0@(comm)
  597. +     .word    0x4a68,0xfff0,0x6bfa    | test
  598. + |    fmulx    fp3,fp6
  599. +     movew    #0x0f23,a0@(comm)
  600. +     .word    0x4a68,0xfff0,0x6bfa    | test
  601. + |    faddx    fp6,fp5
  602. +     movew    #0x1aa2,a0@(comm)
  603. +     .word    0x4a68,0xfff0,0x6bfa    | test
  604. + 
  605. + |    fmulx    fp2,fp4
  606. +     movew    #0x0a23,a0@(comm)
  607. +     .word    0x4a68,0xfff0,0x6bfa    | test
  608. + |    fmulx    fp3,fp0
  609. +     movew    #0x0c23,a0@(comm)
  610. +     .word    0x4a68,0xfff0,0x6bfa    | test
  611. + |    fmulx    fp1,fp2
  612. +     movew    #0x0523,a0@(comm)
  613. +     .word    0x4a68,0xfff0,0x6bfa    | test
  614. + |    fmulx    fp1,fp3
  615. +     movew    #0x05a3,a0@(comm)
  616. +     .word    0x4a68,0xfff0,0x6bfa    | test
  617. + 
  618. + |    faddx    fp3,fp4
  619. +     movew    #0x0e22,a0@(comm)
  620. +     .word    0x4a68,0xfff0,0x6bfa    | test
  621. + |    fdivx    fp5,fp4
  622. +     movew    #0x1620,a0@(comm)
  623. +     .word    0x4a68,0xfff0,0x6bfa    | test
  624. + |    fsubx    fp0,fp2
  625. +     movew    #0x0128,a0@(comm)
  626. +     .word    0x4a68,0xfff0,0x6bfa    | test
  627. + |    fdivx    fp5,fp2
  628. +     movew    #0x1520,a0@(comm)
  629. +     .word    0x4a68,0xfff0,0x6bfa    | test
  630. + 
  631. + 
  632. + |    fmoved    fp4,a1@
  633. +     movew    #0x7600,a0@(comm)        | 
  634. +     .long    0x0c688900, 0xfff067f8
  635. +     movel    a0@,a1@
  636. +     movel    a0@,a1@(4)
  637. + 
  638. + |    fmoved    fp2,a1@(8)
  639. +     movew    #0x7500,a0@(comm)        | 
  640. +     .long    0x0c688900, 0xfff067f8
  641. +     movel    a0@,a1@(8)
  642. +     movel    a0@,a1@(12)
  643. + ");    /* end asm    */
  644. + #endif    sfp004
  645. + 
  646. + #if defined (__M68881__) || defined (sfp004)
  647. + # ifdef ERROR_CHECK    /* no error checking for now    */
  648. + 
  649. + __asm("    
  650. +     pea    _funcname
  651. +     jmp    c_err_check
  652. + ");    /* end asm    */
  653. + 
  654. + # else  ERROR_CHECK
  655. + 
  656. + __asm("rts");
  657. + 
  658. + # endif ERROR_CHECK
  659. + #endif defined (__M68881__) || defined (sfp004)
  660. *** 1.3    1992/02/03 20:19:23
  661. --- cmult.c    1992/03/06 19:38:00
  662. ***************
  663. *** 62,70 ****
  664.    *
  665.    */
  666.   
  667. ! #if defined (__M68881__) && !defined (_M68881)
  668. ! /*# define _M68881*/
  669. ! #endif
  670.   
  671.   #include <stdio.h>
  672.   #include <math.h>
  673. --- 62,68 ----
  674.    *
  675.    */
  676.   
  677. ! #if !defined (__M68881__) && !defined (sfp004)
  678.   
  679.   #include <stdio.h>
  680.   #include <math.h>
  681. ***************
  682. *** 81,83 ****
  683. --- 79,203 ----
  684.   
  685.       return (result);
  686.   }
  687. + #endif !defined (__M68881__) && !defined (sfp004)
  688. + #ifdef    __M68881__
  689. + __asm("
  690. + .text
  691. + .even
  692. + _funcname:
  693. +     .ascii    \"cmult\\0\"
  694. +     .even
  695. + 
  696. + .globl    _cmult
  697. + _cmult:
  698. +     fmoved    sp@(4),fp0
  699. +     fmoved    sp@(12),fp1
  700. +     fmoved    sp@(20),fp2
  701. +     fmoved    sp@(28),fp3
  702. +     fmovex    fp0,fp4
  703. +     movel    a1,d0        | pointer to result
  704. + 
  705. +     fmulx    fp2,fp4
  706. +     fmulx    fp3,fp0
  707. +     fmulx    fp1,fp2
  708. +     fmulx    fp1,fp3
  709. + 
  710. +     fsubx    fp3,fp4
  711. +     faddx    fp0,fp2
  712. + 
  713. +     fmoved    fp4,a1@
  714. +     fmoved    fp2,a1@(8)
  715. + ");    /* end asm    */
  716. + #endif    __M68881__
  717. + 
  718. + #ifdef    sfp004
  719. + __asm("
  720. + 
  721. + comm =     -6
  722. + resp =    -16
  723. + zahl =      0
  724. + 
  725. + .even
  726. + .text
  727. + .even
  728. + _funcname:
  729. +     .ascii    \"cmult\\0\"
  730. +     .even
  731. + .text
  732. + .even
  733. + .globl    _cmult
  734. + _cmult:
  735. + 
  736. +     lea    0xfffa50,a0
  737. + 
  738. +     movew    #0x5400,a0@(comm)    | z1.real -> fp0
  739. +     movel    a1,d0        | pointer to result
  740. +     .long    0x0c688900, 0xfff067f8
  741. +     movel    a7@(4),a0@        | load arg_hi
  742. +     movel    a7@(8),a0@        | load arg_low
  743. + 
  744. +     movew    #0x5480,a0@(comm)    | z1.imag -> fp1
  745. +     .long    0x0c688900, 0xfff067f8
  746. +     movel    a7@(12),a0@        | load arg_hi
  747. +     movel    a7@(16),a0@        | load arg_low
  748. + 
  749. +     movew    #0x5500,a0@(comm)    | z2.real -> fp2
  750. +     .long    0x0c688900, 0xfff067f8
  751. +     movel    a7@(20),a0@        | load arg_hi
  752. +     movel    a7@(24),a0@        | load arg_low
  753. + 
  754. +     movew    #0x5580,a0@(comm)    | z2.imag -> fp3
  755. +     .long    0x0c688900, 0xfff067f8
  756. +     movel    a7@(28),a0@        | load arg_hi
  757. +     movel    a7@(32),a0@        | load arg_low
  758. + 
  759. +     movew    #0x0200,a0@(comm)    | copy fp0 to fp4
  760. +     .word    0x4a68,0xfff0,0x6bfa    | test
  761. + 
  762. + |    fmulx    fp2,fp4
  763. +     movew    #0x0a23,a0@(comm)
  764. +     .word    0x4a68,0xfff0,0x6bfa    | test
  765. + |    fmulx    fp3,fp0
  766. +     movew    #0x0c23,a0@(comm)
  767. +     .word    0x4a68,0xfff0,0x6bfa    | test
  768. + |    fmulx    fp1,fp2
  769. +     movew    #0x0523,a0@(comm)
  770. +     .word    0x4a68,0xfff0,0x6bfa    | test
  771. + |    fmulx    fp1,fp3
  772. +     movew    #0x05a3,a0@(comm)
  773. +     .word    0x4a68,0xfff0,0x6bfa    | test
  774. + |    fsubx    fp3,fp4
  775. +     movew    #0x0e28,a0@(comm)
  776. +     .word    0x4a68,0xfff0,0x6bfa    | test
  777. + |    faddx    fp0,fp2
  778. +     movew    #0x0122,a0@(comm)
  779. +     .word    0x4a68,0xfff0,0x6bfa    | test
  780. + 
  781. + |    fmoved    fp4,a1@
  782. +     movew    #0x7600,a0@(comm)        | 
  783. +     .long    0x0c688900, 0xfff067f8
  784. +     movel    a0@,a1@
  785. +     movel    a0@,a1@(4)
  786. + 
  787. + |    fmoved    fp2,d1
  788. +     movew    #0x7500,a0@(comm)        | 
  789. +     .long    0x0c688900, 0xfff067f8
  790. +     movel    a0@,a1@(8)
  791. +     movel    a0@,a1@(12)
  792. + ");    /* end asm    */
  793. + #endif    sfp004
  794. + 
  795. + #if defined (__M68881__) || defined (sfp004)
  796. + # ifdef ERROR_CHECK    /* no error checking for now    */
  797. + 
  798. + __asm("    
  799. +     pea    _funcname
  800. +     jmp    c_err_check
  801. + ");    /* end asm    */
  802. + 
  803. + # else  ERROR_CHECK
  804. + 
  805. + __asm("rts");
  806. + 
  807. + # endif ERROR_CHECK
  808. + #endif defined (__M68881__) || defined (sfp004)
  809. *** 1.4    1992/02/03 20:19:23
  810. --- csin.c    1992/03/06 19:38:01
  811. ***************
  812. *** 61,69 ****
  813.    *
  814.    */
  815.   
  816. ! #if defined (__M68881__) && !defined (_M68881)
  817. ! /*# define _M68881*/
  818. ! #endif
  819.   
  820.   #include <stdio.h>
  821.   #include <math.h>
  822. --- 61,67 ----
  823.    *
  824.    */
  825.   
  826. ! #if !defined (__M68881__) && !defined (sfp004)
  827.   
  828.   #include <stdio.h>
  829.   #include <math.h>
  830. ***************
  831. *** 78,80 ****
  832. --- 76,190 ----
  833.       result.imag = cos (z.real) * sinh (z.imag);
  834.       return (result);
  835.   }
  836. + #endif !defined (__M68881__) && !defined (sfp004)
  837. + 
  838. + #ifdef    __M68881__
  839. + __asm("
  840. + .text
  841. + _funcname:
  842. +     .ascii    \"csin\\0\"
  843. +     .even
  844. +     .globl _csin
  845. + _csin:
  846. +     movel    a1,d0        | save a1 as return value
  847. +     fmoved    sp@(4),fp0    | z.real
  848. +     fcosx    fp0,fp1        | cos(z.real)
  849. +     fsinx    fp0,fp0        | sin(z.real)
  850. + 
  851. +     fmoved    sp@(12),fp3    | z.imag
  852. +     fcoshx    fp3,fp2        | cosh(z.imag)
  853. +     fsinhx    fp3,fp3        | sinh(z.imag)
  854. +     
  855. +     fmulx    fp2,fp0        | result.real
  856. +     fmoved    fp0,a1@        |
  857. +     fmulx    fp3,fp1        | result.imag
  858. +     fmoved    fp1,a1@(8)    |
  859. + ");    /* end asm    */
  860. + #endif    __M68881__
  861. + 
  862. + #ifdef    sfp004
  863. + __asm("
  864. + | double precision floating point stuff for Atari-gcc using the SFP004
  865. + | developed with gas
  866. + |
  867. + | double precision complex sin
  868. + |
  869. + | M. Ritzert (mjr at dmzrzu71)
  870. + |
  871. + | 12.10.1990
  872. + |
  873. + | addresses of the 68881 data port. This choice is fastest when much data is
  874. + | transferred between the two processors.
  875. + 
  876. + comm =     -6
  877. + resp =    -16
  878. + zahl =      0
  879. + 
  880. + | waiting loop ...
  881. + |
  882. + | wait:
  883. + | ww:    cmpiw    #0x8900,a1@(resp)
  884. + |     beq    ww
  885. + | is coded directly by
  886. + |    .long    0x0c688900, 0xfff067f8
  887. + | and
  888. + | www:    tst.b    a1@(resp)
  889. + |    bmi.b    www
  890. + | is coded by
  891. + |    .word    0x4a68,0xfff0,0x6bfa        | test
  892. + 
  893. +     .text; .even
  894. +     .globl _csin
  895. + _csin:
  896. +     movel    a1,d0                | save a1 as return value
  897. +     lea    0xfffa50,a0            | fpu address
  898. + 
  899. +     movew    #0x5431,a0@(comm)        | sincos: sin -> fp0 
  900. +     .long    0x0c688900, 0xfff067f8        |      cos -> fp1
  901. +     movel    sp@(4), a0@            | load z.real
  902. +     movel    sp@(8), a0@            | load z.real
  903. + 
  904. + |    fmoved    sp@(12),fp3            | z.imag to fp3
  905. +     movew    #0x5580,a0@(comm)
  906. +     .long    0x0c688900, 0xfff067f8
  907. +     movel    sp@(12),a0@
  908. +     movel    sp@(16),a0@
  909. + 
  910. + |    fcoshx    fp3,fp2                | cosh(z.imag)
  911. +     movew    #0x0d19,a0@(comm)
  912. +     .word    0x4a68,0xfff0,0x6bfa        | test
  913. + |    fsinhx    fp3,fp3                | sinh(z.imag)
  914. +     movew    #0x0d82,a0@(comm)
  915. +     .word    0x4a68,0xfff0,0x6bfa        | test
  916. + |    fmulx    fp2,fp0                | result.real
  917. +     movew    #0x0823,a0@(comm)
  918. +     .word    0x4a68,0xfff0,0x6bfa        | test
  919. + |    fmoved    fp0,a1@                |
  920. +     movew    #0x7400,a0@(comm)        | 
  921. +     .long    0x0c688900, 0xfff067f8
  922. +     movel    a0@,a1@
  923. +     movel    a0@,a1@(4)
  924. + |    fmulx    fp3,fp1                | result.imag
  925. +     movew    #0x0ca3,a0@(comm)
  926. +     .word    0x4a68,0xfff0,0x6bfa        | test
  927. + |    fmoves    fp1,d1                |
  928. +     movew    #0x7480,a0@(comm)        |
  929. +     .long    0x0c688900, 0xfff067f8
  930. +     movel    a0@,a1@(8)
  931. +     movel    a0@,a1@(12)
  932. + ");    /* end asm    */
  933. + #endif    sfp004
  934. + 
  935. + 
  936. + #if defined (__M68881__) || defined (sfp004)
  937. + # ifdef ERROR_CHECK    /* no error checking for now    */
  938. + __asm("
  939. +     pea    _funcname
  940. +     jmp    c_err_check
  941. + ");    /* end asm    */
  942. + # else  ERROR_CHECK
  943. + 
  944. + __asm("rts");
  945. + 
  946. + # endif ERROR_CHECK
  947. + #endif defined (__M68881__) || defined (sfp004)
  948. *** 1.3    1992/02/03 20:19:23
  949. --- csubt.c    1992/03/06 19:38:02
  950. ***************
  951. *** 124,131 ****
  952.   
  953.       movew    #0x5428,a0@(comm)    | fp0 -= z2.real
  954.       .long    0x0c688900, 0xfff067f8
  955. !     movel    a7@(4),a0@        | load arg_hi
  956. !     movel    a7@(8),a0@        | load arg_low
  957.   
  958.       movew    #0x5480,a0@(comm)    | z1.imag -> fp1
  959.       .long    0x0c688900, 0xfff067f8
  960. --- 124,131 ----
  961.   
  962.       movew    #0x5428,a0@(comm)    | fp0 -= z2.real
  963.       .long    0x0c688900, 0xfff067f8
  964. !     movel    a7@(20),a0@        | load arg_hi
  965. !     movel    a7@(24),a0@        | load arg_low
  966.   
  967.       movew    #0x5480,a0@(comm)    | z1.imag -> fp1
  968.       .long    0x0c688900, 0xfff067f8
  969. ***************
  970. *** 133,142 ****
  971.       movel    a7@(16),a0@        | load arg_low
  972.   
  973.       movew    #0x54a8,a0@(comm)    | fp1 -= z2.imag
  974. !     movel    a1,d0             pointer to result
  975.       .long    0x0c688900, 0xfff067f8
  976. !     movel    a7@(4),a0@        | load arg_hi
  977. !     movel    a7@(8),a0@        | load arg_low
  978.   
  979.   |    movew    #%0111 0101 0000 0000,a0@(comm)    | fetch fp0
  980.       movew    #0x7400,a0@(comm)        | 
  981. --- 133,142 ----
  982.       movel    a7@(16),a0@        | load arg_low
  983.   
  984.       movew    #0x54a8,a0@(comm)    | fp1 -= z2.imag
  985. !     movel    a1,d0            | pointer to result
  986.       .long    0x0c688900, 0xfff067f8
  987. !     movel    a7@(28),a0@        | load arg_hi
  988. !     movel    a7@(32),a0@        | load arg_low
  989.   
  990.   |    movew    #%0111 0101 0000 0000,a0@(comm)    | fetch fp0
  991.       movew    #0x7400,a0@(comm)        | 
  992.